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Abstract 

This paper presents an a priori error analysis of the /ip- version of the boundary element 
method for the electric field integral equation on a piecewise plane (open or closed) Lipschitz 
surface. We use H(div)-conforming discretisations with Raviart-Thomas elements on a se- 
quence of quasi- uniform meshes of triangles and/or parallelograms. Assuming the regularity 
of the solution to the electric field integral equation in terms of Sobolev spaces of tangential 
vector fields, we prove an a priori error estimate of the method in the energy norm. This 
estimate proves the expected rate of convergence with respect to the mesh parameter h and 
the polynomial degree p. 
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1 Introduction 

With this paper we continue the analysis of high-order boundary element methods (BEM) for the 
electric field integral equation (EFIE) started in [31[6]. Our BEM is based on discretisations of 
the variational formulation of the EFIE (called Rumsey's principle) with an H(div)-conforming 
family of boundary elements. This approach is referred to as the natural boundary element 
method for the EFIE. In [4] we analysed the natural p-BEM for the EFIE on a plane open 
surface with polygonal boundary. We proved convergence of the p- version with Raviart-Thomas 
(RT) parallelogram elements and derived an a priori error estimate which takes into account the 
strong singular behaviour of the solution at edges and corners of the surface. In our previous 
paper [6j we considered the EFIE on a piecewise plane (open or closed) Lipschitz surface P and 
proved quasi-optimal convergence of the natural /ip-BEM with quasi-uniform meshes of triangles 
and quadrilaterals. In the present note we perform an a priori error analysis of that method on 



*Supported by EPSRC under grant no. EP/E058094/1. 

^Department of Mathematical Sciences, Brunei University, Uxbridge, West London UBS 3PH, UK. Email: 
albespalovOyahoo . com 

^Facultad de Matematicas, Pontificia Universidad Catolica de Chile, Avenida Vicuria Mackenna 4860, Santiago, 
Chile. Email: nheuer@mat . puc . cl 



1 



affine meshes under the assumption that the regularity of the exact solution is given in Sobolev 
spaces of tangential vector fields on F. As the main result we prove an a priori error estimate in 
the energy norm (Theorem 12. 2p . The estimate appears to be optimal with respect to the mesh 
size h and the polynomial degree p as the convergence rates in both h and are r + 1/2 for p 
large enough. This corresponds to the expected rate which is the Sobolev regularity order r of 
the exact solution minus the Sobolev order —1/2 of the energy norm. 

While in the /i- version the degrees of approximating polynomials are fixed (usually at a low 
level) and convergence is achieved by refining the mesh, the p-version keeps the mesh fixed 
and improves approximations by increasing polynomial degrees. The /ip- version combines both 
mesh refinement and increase of polynomial degrees. For boundary integral equations governing 
the Laplace equation optimal /ip-convergence rates for singular problems (and quasi-uniform 
meshes) are proved in [5l [3] . The analysis of optimal /ip-BEM convergence rates for the EFIE 
with singular solutions is an open problem and under investigation. In this paper we deal with 
the case of solutions with Sobolev regularity. We also note that the /ip-BEM with geometrically 
graded meshes (yielding an exponential rate of convergence) has been studied in [23], again 
for the Laplace equation and for hypersingular and weakly singular integral operators. For the 
EFIE its analysis is an open problem. 

An a priori error analysis of the natural /i-BEM for the EFIE was performed in [25] for 
polyhedral surfaces and in [11] for open Lipschitz surfaces (see also jl5j for a survey of results 
and techniques). In particular, an optimal /i-convergence rate of the method for given Sobolev 
regularity of the solution and given polynomial degree has been proved in [11] (being proved 
on open surfaces, this result extends to polyhedral surfaces as well). If the BEM for the EFIE 
converges quasi-optimally, then a priori error analysis reduces to an approximation problem 
within the energy space, which is either H~^/^(divr, F) or Hq ^'''^(divr, F) depending on whether 
the surface F is closed or open. The main tool in the approximation analysis is an appropriate 
H(div)-conforming interpolation operator. Whereas the standard H(div)-conforming RT inter- 
polation operator works well for the /i- version (cf. [llj). it does not provide an optimal result for 
high order methods. The main reason for this is the lack of stability (with respect to polynomial 
degrees) of this operator for low-regular vector fields, which always appear when dealing with 
the EFIE on non-smooth surfaces. Furthermore, existing techniques to prove p-estimates for the 
error of RT interpolation work only on quadrilateral elements, and their extension to triangular 
elements does not seem feasible. 

An alternative to the classical RT interpolation operator is a corresponding projection based 
interpolation operator. Such operators were first introduced in [19] to analyse high-order finite 
element approximations with H(curl)-conforming edge elements for Maxwell's equations in two 
dimensions. However, applying a simple rotation argument the results of \T9\ can be formulated 
in the H(div)-conforming setting, which is intrinsic to natural boundary element discretisations 
of the EFIE (see also [7]). The projection based interpolation operators are stable with respect 
to polynomial degrees, they work equally well on both triangular and quadrilateral elements 
and also for low-regular fields. That is why these operators have become an efficient tool in the 
analysis of high-order methods (see [Oj [U [24] for the finite element methods and [4t |6| for the 
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BEM). The a priori error analysis in this paper relies on such an operator as well. In particular, 
we demonstrate that employing the H(div)-conforming projection based interpolation operator 
(rather than the classical RT interpolation operator) one obtains an optimal error estimate for 
the /ip-BEM with quasi- uniform meshes. 

The paper is organised as follows. In the next section we formulate the EFIE (in a variational 
form) and define the /ip- version of the BEM with quasi- uniform meshes. We also formulate the 
main result (Theorem [22]), which states an a priori error estimate of the approximation method. 
Section [3] gives necessary preliminaries: first, in ^S.ll we introduce the needed notation and recall 
definitions of Sobolev spaces of tangential vector fields; then, in §3.21 we sketch the definition 
of the H(div)-conforming projection based interpolation operator on the reference element and 
prove a new property of this operator related to approximations of normal traces on the element's 
edges (Lemma I3.3p . In Section U] we study approximating properties of the discrete (boundary 
element) space X/ip in the energy space X of the EFIE, and prove that the orthogonal projection 
onto X/ip with respect to the norm in X satisfies an optimal error estimate in both h and p. 

Throughout the paper, C denotes a generic positive constant which is independent of h, p 
and involved functions. 

2 Formulation of the problem and the main result 

We consider the EFIE for which we have proved quasi-optimal convergence of the /ip-BEM with 
quasi- uniform meshes [6j. In this paper we provide an a priori error estimate. To this end let 
us recall the model problem, its /ip-discretisation and the involved spaces. 

Let r denote a piecewise plane (open or closed) Lipschitz surface in W^. In the case of an 
open surface we additionally assume that F is orientable. Let us introduce Rumsey's formulation 
of the electric field integral equation on F. For a given wave number k > and a scalar function 
V (resp., tangential vector field v) we define the single layer operator \I'fc (resp., ^k) by 

1 f Qik\x-y\ 

<l!kv{x) = — / v{y)- dSy, X £ R\r 

4vr 7r \x - y\ 

resp., *fev(x) = — / ^r{y)- -dSy, x G MV . 

Let Lj2(F) be the space of two-dimensional, tangential, square integrable vector fields on F. 
By Vr (resp., divr) we denote the surface gradient (resp., surface divergence) acting on scalar 
functions (resp., tangential vector fields) on F. We will need the following space: 

X = H-i/2(divr,F) := {u G H^^/^(F); divr u G R-^/^F)} 

if F is a closed surface, and 

X = Ho^/^(divr,F) := {u G H^^^^(F); divr u G H-^/^{T) and 

(u, Vrv) + (divr u, w) = for all v G C7°°(F)} 
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if r is an open surface. In the latter definition the brackets (•, •) denote dualities associated with 
H|'|''^(r) and H^/'^{T), respectively. For definitions of the space C°°{T) and the Sobolev spaces 
on r we refer to ^3.11 below. Throughout, we use boldface symbols for vector fields. The spaces 
(or sets) of vector fields are also denoted in boldface (e.g., H*(r) = {H'^{T))^). 

Let X' be the dual space of X (with L^(r) as pivot space). Now, for a given tangential vector 
field f G X' (f represents the excitation by an incident wave), Rumsey's formulation reads as: 
find a complex tangential field u G X such that 

a(u,v) := (7tr(^fcdivru),divrv) - A:2(,r,(*^u),v) = (f,v) Vv G X. (2.1) 

Here jtr is the standard trace operator, and vr^ denotes the tangential components trace mapping 
(see ^3.11 for the definition) . To ensure the uniqueness of the solution to (j2.ip in the case of the 
closed surface T we always assume that k"^ is not an electrical eigenvalue of the interior problem. 

For the approximate solution of (j2.ip we apply the /ip- version of the BEM based on Galerkin 
discretisations with Raviart-Thomas spaces on quasi-uniform meshes. In what follows, h > 
and p > 1 will always specify the mesh parameter and a polynomial degree, respectively. For 
any C we will denote pQ = sup{diam(i3); B is a ball in 0,}. Furthermore, throughout the 
paper, K is either the equilateral reference triangle T = {x2 > 0, X2 < xiVs, X2 < (1 — xi)\/^} 
or the reference square Q = (0, 1)^. A generic side of K will be denoted by i. 

Let T = {Afi} be a family of meshes = {Tj; j = 1, . . . , J} on F, where the elements 
Tj are open triangles or parallelograms such that f = U^^^Fj, and the intersection of any two 
elements Tj, F^ {j ^ k) is either a common vertex, an entire side, or empty. 

We denote hj = diam(Fj) for any Tj G A/j. The elements are assumed to be shape regular, 

i.e., there exists a positive constant C independent oi h = max hj such that for any Tj G A/j 

j 

and arbitrary A/j G T there holds hj < C pr^ ■ Furthermore, any element Tj is the image of the 
corresponding reference element K under an affine mapping Tj, more precisely 

f,- = T,{K), X = r,(^), X = {xi,x2) G f,-, ^ = (6,6) e k. 

The Jacobian matrix of Tj is denoted by DTj and its determinant Jj := det{DTj) satisfies the 
relation \Jj\ — h'j. 

We consider a family T of quasi-uniform meshes A/j on F in the sense that there exists a 
positive constant C independent of h such that for any Tj G A/j and arbitrary A/j G T there 
holds h < C hj. 

The mapping Tj introduced above is used to associate the scalar function u defined on the 
real element Tj with the function u defined on the reference element K: 

u = uo T^^ on Tj and u = uoTj on K. 

Any vector-valued function v defined on K is transformed to the function v on Fj by using the 
Piola transformation: 

V = Mj{w) = j-DTjV o Tr\ v = Mj\v) = JjDTf^v o T,-. (2.2) 
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Let us introduce the needed polynomial sets. By Vp{I) we denote the set of polynomials 
of degree < p on an interval / C M, and V^^I) denotes the subset of T'p{I) which consists of 
polynomials vanishing at the end points of /. In particular, these two sets will be used for an 
edge i C dK. 

Further, 'Pp(r) denotes the set of polynomials on T of total degree < p, and 'Ppi^p^iQ) is the 
set of polynomials on Q of degree < pi in ,^1 and degree < p2 in ^2- For pi = p2 = p we denote 
'PpiQ) = 'Pp,piQ), '^ill ^se the unified notation T'p{K), which refers to Vp{T) \i K = T 

and to Vp{Q) if K = Q. The corresponding set of polynomial (scalar) bubble functions on K is 
denoted by V^{K). 

Let us denote by Tf^{K) the RT -space of order p > 1 on the reference element K (see, e.g., 

mm), i-e-, 



Vf'{K) = {Vp-,{K)r®^Vp^,{K) 



{Vl_,{T)f®$.Vl^,{T) ii K = T, 

Kp-i(Q)''K~iAQ) '^^ k = q. 



The subset of V^"^ {K) which consists of vector- valued polynomials with vanishing normal trace 
on the boundary dK (vector bubble- functions) will be denoted by ^^^'^{K). 
Then using transformations ()2.2I) . we set 

l^hp ■■= {v G XO; MJH^r\r^) G Vf^{K), j = 1, . . . , J}, (2.3) 

where the space X'^ C X is defined in ^3.11 (X'^ = H(divr, T) if F is closed and X'^ = Ho(divr, F) 
if F is an open surface). We will denote by = N{h,p) the dimension of the discrete space 
X/jp. One has ~ for fixed p and N c:^ p"^ for fixed h. 

The hp-version of the Galerkin BEM for the EFIE reads as: Find u^p £ X/^p such that 

a(u;,p,v) = (f,v) VvGX;,p. (2.4) 



First, let us formulate the result which states the unique solvability of (j2.4p and quasi-optimal 
convergence of the hp-veision of the BEM for the EFIE. 

Theorem 2.1 [6, Theorem 2.1] There exists Nq > 1 such that for any f G X' and for arbitrary 
mesh-degree combination satisfying N{h,p) > A'^o the discrete problem ()2.4p is uniquely solvable 
and the hp-version of the Galerkin BEM generated by RT-elements converges quasi- optimally, 
i.e., 

||u - u/ip||x < Cinf{||u - v||x; v G X/,p}. (2.5) 

Here, u G X is the solution of (j2.ip . u/jp G X/jp is the solution of (|2.4p . || • ||x denotes the norm 
in X, and C > is a constant independent of h and p. 

The following theorem is the main result of this paper (its formulation involves the space X^ 
and the norm 11 • llx'- which are defined in ^3.ip . 
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Theorem 2.2 Let u G X and Uhp £ X/jp be the solutions of ()2.ip and ()2.4p . respectively. Then 
there exists a real number r (r>OifTis closed, and r G (— 1,0) i/F is an open surface) such 
that u G X^ and the following a priori error estimate holds 

||u-U;,p||x < C/ll/2+min{r,p}^-{r+l/2) ||^||^^ (2.6) 

with a positive constant C independent of h and p. 

Proof. The assertion regarding the regularity of the solution u to the EFIE is a direct con- 
sequence of the regularity results in [XT', Section 4.4] (see also [H Appendix A]). Due to the 
quasi-optimal convergence (j2.5p of the /ip-BEM, the error estimate in ()2.6p then immediately 
follows from the approximation result of Theorem 14.11 below. □ 

Remark 2.1 If T is a closed Lipschitz polyhedral surface (T = d^) and the EFIE represents 
a boundary value problem for the time-harmonic Maxwell's equations in the exterior domain 
then one can be more specific about the regularity of the solution u of (I2.ip . In fact, the 
presence of re-entrant edges is inevitable for the exterior of any polyhedral domain. Therefore, 
using the results for Maxwell singularities (see ^\), we conclude that for a closed polyhedral 
surface T there exists r G (0, ^) such that u G X'', and u 'S}/'^ . 

Remark 2.2 The convergence rate of the hp-BEM for the EFIE is limited by the low regularity 
of the solution to the corresponding Maxwell's equations, especially for problems with screens 
which represent the least regular case. To obtain sharp a priori error estimates, a refined ap- 
proximation analysis of singularities inherent to the solution of the EFIE is needed (see j^J. 
However, when the regularity of the solution to the EFIE is stated in terms of Sobolev spaces of 
tangential vector fields on T, the result of Theorem \2.2\ is optimal with respect to both h and p. 

3 Preliminaries 

3.1 Functional spaces, norms, and inner products 

In our previous paper [6J we recalled definitions of the full range of Sobolev spaces necessary 
for the convergence analysis of the BEM for the EFIE (see §3.1 therein). This included Sobolev 
spaces on a Lipschitz domain C M" and Sobolev spaces of scalar functions and tangential 
vector fields on a piecewise smooth (open or closed) Lipschitz surface F C M^. We have also 
defined basic differential operators on F. In the present paper we will use the same notation as in 
[6] for all differential operators, Sobolev spaces and their norms. For convenience of the reader, 
let us repeat some essential definitions, in particular, those for Sobolev spaces of tangential vector 
fields. Furthermore, we will introduce some more spaces and norms, which are indispensable for 
the error analysis of the /ip-BEM. 

Let be a Lipschitz domain in M". We use a traditional notation for the Sobolev spaces 
H'{n) (s > -1), H^{n) {s G (0,1]), and H'{n) (s G [-1,1]) with their standard norms 
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(cf., [26j). In particular, for s £ (0,1), the spaces H''{Q) are defined by using the real K- 
method of interpolation, and for s £ [—1,0), the spaces H^{^}), H^{^}) and their norms are 
defined by duality with L'^{Q) = H^{^) = H^{^) as pivot space. The norm and inner product 
in L^(J7) will be denoted as || • ||o,r2 and (•, •)o,q, respectively. 

It is known that the standard norms || • and || • are not scalable for s G (0, 1] 

under affine transformations of 0, onto the reference domain (element). However, we will need a 
scalable norm in the space on a generic interval £ C M^. Scalable families of norms in H'^ 
and in H^^ for s £ [—1,1] were introduced in [20]. Following [20] we define 



\Hlie) = (meas(^)) ||n||o,£ + 
Then the norm || • H/^-i^^) is defined by duality: 

Let £ be the image of the reference interval i under an affine transformation M, i.e., i = M{1), 
and let meas(£) ~ h. Denote f = f o M. Then the norms || • and || • Hj^^-i^^^ are scalable 

(see pni Lemma 3.1]), i.e., 

for any / G H^{i) and / G H~^{£), respectively, and both equivalences are uniform for h> 0. 
An important fact related to the norm || • \\g-i is that it enjoys the localisation property, 

h 

provided that the function has zero average on each sub-domain (see [201 Lemma 3.2]). In 
particular, we will need the following result. 



Lemma 3.1 Assume that the boundary dQ, of the polygonal domain J7 C is partitioned into 
N segments ij (j = 1,...,N). Then, for all f £ H'^{dn) with f\i. £ H'^iij) and f da = 
(j = 1, . . . , N), there holds 



N 



with a positive constant C independent of f and N. 

Now, let r be a piecewise smooth (open or closed) Lipschitz surface in M^. We will assume 
that r has plane faces T^^^ (i = 1, . . . ,X; without loss of generality it is assumed that X > 1) 
and straight edges Cij = r« nr(j) / (i / j). if r is a closed surface, we will denote by $7 
the Lipschitz polyhedron bounded by T, i.e., F = d^l. If F is an open surface, we additionally 
assume that F is orientable. In this case, we first introduce a piecewise plane closed Lipschitz 
surface F which contains F, and then denote by O the Lipschitz polyhedron bounded by F, i.e., 
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f = d^l. For each face F^*^ C F there exists a constant unit normal vector i/j, which is an outer 
normal vector to 0,. These vectors are then blended into a unit normal vector v defined almost 
everywhere on F. For each pair of indices i,j = 1, . . . ,1 such that F^*-* n F^-'-' = eij we consider 

unit vectors Tij, , and Tj such that Tij\\eij, t- = Tij x Vi, and Tj = Tij x Vj. Since each 

F^*) can be identified with a bounded subset in M^, the pair (T^^Tij) is an orthonormal basis 
of the plane generated by F^*). 

Let F be a closed surface. Then F = dfl is locally the graph of a Lipschitz function. 
Since the Sobolev spaces -ff* for |s| < 1 are invariant under Lipschitz (i.e., C^'^) coordinate 
transformations, the spaces H^{T) with |s| < 1 are defined in the usual way via a partition of 
unity subordinate to a finite family of local coordinate patches (see [27j). Due to this definition, 
the properties of Sobolev spaces on Lipschitz domains in M" carry over to Sobolev spaces on 
Lipschitz surfaces. If F is an open surface, then the Sobolev spaces H^{T), H^{T) for |s| < 1 
and Hq{T) for < s < 1 are constructed in terms of the Sobolev spaces H^(T) on a closed 
Lipschitz surface f D F (see [27]). Note that the spaces i7*(F(*)) and ^^(F^*)) on each face F^*) 
are well-defined for any s > — 1. 

We will denote by jtr the standard trace operator, jtriu) = u\r, u E C°°(0). For s S (0, 1) 
(resp., s > 1), 7tr has a unique extension to a continuous operator H'^~^^/'^{i}) —>■ H^{T) (resp., 
H'+^/^{Q) H^{r)), see [Ml H]. We wiU use the notation C°^{T) = -/tr{C°^{n)). 

Using the introduced Sobolev spaces of scalar functions, we define: 

= {H'{n)f, H^(F») = (i7"(FW))^ ff(FW) = (i?"(FW))2 (1 < i < X) 
for s > —1, and 

ll'{T) = {H'{T)f for sG [-1,1]. 

The norms and inner products in all these spaces are defined component- wise and usual conven- 
tions H0(J7) = L'^in), H0(F) = L2(F), H0(F«) = H0(F») = L2(F») hold. 

Now let us introduce the Sobolev spaces of tangential vector fields defined on F (see [121 13 
[14]). We start with the space 

L?(F) := {u G L2(F); u • = on F}, 

which will be identified with the space of two-dimensional, tangential, square integrable vector 
fields. The norm and inner product in this space will be denoted by || • ||o,r and (•, •)o,ri respec- 
tively. The similarity of this notation with the one for scalar functions should not lead to any 
confusion. Then we define (hereafter, Uj denotes the restriction of u to the face F^*)): 

H^_(F) := {u G L2(F); Ui G H^(F«), 1 < i < J}, s > 0, 



lullH-.cr) := I E l|uillH=(r{0) 

\i=l 



Let 7tr be the trace operator (now acting on vector fields), 7tr(u) = u|r, 7tr : H*+-'^/^(r2) — > 
H'^(F) for s G (0,1), and let 7^^"*^ be one of its right inverses. We will use the "tangential 
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components trace" mapping tTt : {C°°(Q))^ L^(r) and the "tangential trace" mapping 7,- : 
(C°°(J1))^ — > L^(r), which are defined as u 1/ x (u x i>')|r and u x i/|r, respectively. We will 
also use the notation tTt (resp., 7,-) for the composite operator vr,- o 7^"^ (resp., jr ° Itr )' which 
acts on traces. Then we define the spaces 

H;/'(r) := 7r.(HV2(r)), nf (F) := 7r(HV2(r)), 

endowed with their operator norms 

l|u|lHi/2fn •= ^ ™^ {ll0llHl/2(r); ^r(0) = u}, 

\Mn^/^(r) '■= ^ {||0llHi/2(r); 7r(<^) = u}. 

^'-> 0eHi/2(r) ^ ^ 

1/2 1/2 

It has been shown in [12] that the space H|| (F) (resp., Hj_ (F)) can be characterised as 

1/2 

the space of tangential vector fields belonging to H _ (F) and satisfying an appropriate "weak 
continuity" condition for the tangential (resp., normal) component across each edge eij of F. 
For s > ^ we set 

H||(F) := {u G H*_(F); Uj • Tij = Uj • Tij at each Cij}, 
Hi(F) := {uGff_(F); u, •rp^ = uj -rf at each 6^^}. 

For any s > ^ the spaces Hy(F) and H^(F) are closed subspaces of H*„(F). Finally, for s G [0, ^) 
we set 

H^(F) = Hi(F) := H^_(F). 

If F is an open surface, then we also need to define subspaces of H|(F) and H^(F) incorpo- 
rating boundary conditions on dT (for tangential and normal components, respectively). In this 
case, for a given function u on F, we will denote by u the extension of u by zero onto a closed 
Lipschitz polyhedral surface F D F. Then we define the spaces 

Hjj(F) := {u G H|(F); a G H|(f)}, s > 0, 
H1(F) := {uGHi(F); uGHi(f)}, s>0, 
which are furnished with the norms 

ll'^llHj|(r) •= ll"llH|{f)' ll"llH^(r) •= ll"llH^(f)' s>0. 

When considering open and closed surfaces at the same time we use the notation Hj|(F), Hj_(F), 
etc. also for closed surfaces by assuming that Hj|(F) = Hj|(F), H^(F) = H^(F), etc. in this 
case. This in particular applies to the following definition of dual spaces. For s G [—1,0), the 
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spaces Hjj(r), H||(r), H^(r), and H5_(r) are defined as the dual spaces of H|['(r), H||"'(r), 

H][''(r), and Hj''(r), respectively (with Lj(r) as pivot space). They are equipped with their 
natural (dual) norms. Moreover, for any s € ( — ^, ^) there holds (cf. [21j ) 

Hf(r) = H|(r) = Hi(r) = Hi(r). 

Using the above spaces of tangential vector fields, one can define basic differential operators 
on r. The tangential gradient, Vr : H^(r) L^(r), and the tangential vector curl, curlp : 
H^{T) — > L^(r), are defined in the usual way by localisation to each face r^*^ The adjoint 
operator of — Vr is the surface divergence denoted by divp (we refer to |12l [T3] for more details 
regarding definitions and properties of differential operators on both closed and open surfaces). 

Now we can introduce the spaces which appear when dealing with the EFIE on T. First, we 

set 

H"(divr,r) := {u G H|(r); divp u G F"(r)}, sG [-1/2,1/2] 

and 

H"_(divr,r) := {u G H"_(r); divp u G -f^'_(r)}, s > 

for r being either a closed or an open surface. Here, the space H'^__ (F) is defined similarly to the 
space H*_(F) in a piecewise fashion: 

H'_{T) := {u G L2(F); -ulr^ G /7^(r»), i = l,... ,1}, s > 0, 

I 

II l|2 II I ||2 

1=1 

If F is an open surface, then we will also use the space 

H"(divr,F) := {u G H|(F); divp u G -H'"(F)}, sG [-1/2,0]. 

The spaces H*(divr,F), H*_(divr,F), and H*(divr,F) are equipped with their graph norms 
II • llH^Cdivr.r), II • llH'>_(divr,r), and || • ||H^(divj,^r)' respectively. If s = 0, then we will drop the su- 
perscript and for open surfaces also the tilde in the above notation, H''(divr, F) = (divp, F) = 
HO(divr,F) =H(divr,F). 

On open surfaces, one needs the spaces incorporating homogeneous boundary conditions for 
the trace of the normal component on dT. By HQ(divr,F) with s G [—^,0] we denote the 
subspace of elements u G H'*(divr,F) such that for all v G C°°(F) there holds 

(u, Vr?;) + (divr u, v) = 0, 

where brackets (•,•) denote the corresponding dualities. Similarly as above, we will drop the 
superscript and the tilde if s = 0. We note that HQ(divr,F) is a closed subspace of H*(divr,F) 
for s G [-^,0]. 
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The notation above related to open surfaces will be used also for two-dimensional domains (in 
particular, for single faces of T and for reference elements). When we need to join the notation 
for open and closed surfaces, we will write 



X = H" 



V2(divr,r), X^ = < 



'H*(divr,r) forsG[-i,0) 

H(divr,r) fors = 0, 
.H^_(divr,r) fors>0 



if r is a closed surface, and 



X = Ho'/'(divr,r 



'Hg(divr,r) forsG[-i,0) 
X" = <^ Ho(divr,r) for s = 0, 

, H*_(divr,r) nHo(divr,r) for s > 



if r is an open surface. In all these cases the norms will be denoted as || • ||x and || • Hx"- 

Let [•, -Je {9 £ [0, 1]) denote the standard interpolation (see [261 [2]). We quote the following 
interpolation result from [11] (see Theorem 4.12 therein). 

Lemma 3.2 There exists an sq S (0, ^] such that for any s S [— ^,,so) there holds [X,X^]fii = 
X(i/2+s)e-i/2_ 

3.2 Interpolation operators 

As it was mentioned in the introduction, our analysis of /ip-approximations essentially relies on 
the properties of the H(div)-conforming projection based interpolation operator. Let us briefly 
sketch the definition of this operator on the reference element K (see [19] for details). In this 
sub-section we use standard differential operators V, curl and div acting on 2D scalar functions 
and vector fields, respectively. 



Given a vector field u G W{K) n H(div, K) with r > 0, the interpolant = IIp^u G 
T'?^{K) is defined as the sum of three terms: 



where fi denotes the outward normal unit vector to dK, and cf)f> are the standard basis functions 
(associated with edges £) for 7-'Y^{K). For any edge i C dK one has 



Ui + + u|. 



Here, ui is the lowest order interpolant 





(3.3) 



e 
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Hence, there exists a scalar function ip, defined on dK, such that 

dip 

— — = (u — ui) • fi, ip = at ah vertices. (3-4) 
da 

Then, for each edge £, the restriction ipli is projected onto the set of polynomials 'Pp^i{£) 

Extending 1^2/ by zero from i onto dK (and keeping its notation) and using the polynomial 
extension from the boundary, we define V'2,p+i ^ "Pp+ii^) such that ip2,p+i\dK = Then we 
set 

"2= ^le^K^i^)^ where u^^ = curl V^p+i- 

The interior interpolant Ug is a vector bubble function that solves the constrained minimization 
problem 

||div(u - (ui + + uf))||o,/f ^ min, 

(u - (ui + + uf ), curl (/<)o,i^ = V,^ G 'Pp+i(K). 

Remark 3.1 The Soholev space H^/"^ and the corresponding norm were defined in [6j using the 
real K-method of interpolation (see [26jj. However, the expression for the ^1/2 

-inner product in 

(j3.5p is based on another (equivalent) norm in H^^^. Without loss of generality, let us assume 
that i = I = (0,1) is the reference interval. Then one has (see, e.g., [19j ) 

UWhi/^I) - l|V^||o,T, 

where T is the reference triangle, (ff denotes the extension of the function cj) £ H^^'^{I) by zero 
onto dT D /, and is the harmonic lift of (ff G H^^^{dT). Then, applying the parallelogram 
law and integrating by parts, we find the expression of the corresponding inner product (cf. [19\): 

\ on / 0,1 \ an / 0,1 
This expression can also be written as 

{(PA)m/2(^i) = {(P:'DMT{i^))^^ = (vMT{<P),i)^^ V(/<, V e ^^/'(/), (3.6) 

where V>Mt '■ H'^^'^{I) H^^^^{I) denotes the Dirichlet-to-Neumann operator associated with 
the triangle T. Exactly this inner product is employed in (13. 5p . 
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For r > the operator 11^'^ : 1¥'{K) n H(div, K) H(div, K) is well defined and bounded, 
with corresponding operator norm being independent of the polynomial degree p (cf. |19l Propo- 
sitions 2]). Moreover, Hp" preserves polynomial vector fields in 7^^'^(K), and there holds the 
following estimate for the interpolation error (see [71 Theorem 5.1]) 



I|U - n^'^ u||H(div,/f) < C'P l|u||H'-(div,i^); (3-7) 

provided that u G H^(div, K) := {u G W{K)- divu G H''{K)} for r > 0. 
We will need the following auxiliary result regarding the operator Hp" ■ 



Lemma 3.3 Let u G W{K) D H(div, K) with r > 0, and let uP = Uf'u G Vp (K). Then /< 
any edge i C dK there holds 

II (u - uP) . n\\H-Hi) < CP''^^ IK" - ■ ^H-^HdK) < Cp-'^^ ||u - u*'||^,(div,i^). (3i 
Proof. Let us fix an edge i C dK. Using the definition of the interpolant we have 



or 



(u - uP) ■ n||j^_i(^) = ||(u - ui) • n - • n||^_i(^) 









9^ 


da 


H 



(3.9) 



Since the derivative with respect to the arc length ^ is a bounded operator mapping L'^{i) to 
H'^ii) (see, e.g., [221 Lemma 3]), we deduce from (|3.9p 



||(u-uP)-n||^_i(,) <C||V|^-^2/||o,£- (3.10) 

Let T be an equilateral triangle having i as one of its edges {if K = T then T = T). Denoting 
(j) '■= — ip2/ and using the expression for the H^^'^{l)-mnei product as in (j3.6p . we can rewrite 
the orthogonality relation in (j3.5p as 



'p/m/He) = \<P^^^fi^p))^, = ^ V0pGPpV(£). (3.11) 

Now, let us consider an auxiliary mixed boundary value problem on T: find ^ G H^{T) such 
that 

-A^> = in r, — = (j) on£, ^> = on dT\£. 
on 

Recalhng that </> G H^/^{£) C H^/^{i) and using the regularity theory for elliptic problems in non- 
smooth domains (see, e.g., [I8l[21]), we conclude that <^ G H'^{T) and ||<I>||j:^2('f) < C* ||0||j:^i/2(£). 
Therefore, applying the trace theorem for a single edge £ C T (see |2H Theorem 1.4.2]) we prove 
that $|f G H^/'^{£) and 

m\H^/2(^e) < C ||^>||^2(f) < C U\\H^/2^gy (3.12) 
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On the other hand, (/> = d^/dn = 'DMf{^\e) and using (j3.1ip we obtain for any c/ip G 'Pp{i) 



Hence, using the fact that the operator VMf is continuous as a mapping Hq{£) L'^i^), we 
find 

m\o,i<C\^U-cl)p\H^^,) V0pGpO(£). (3.13) 

To estimate the semi-norm in (j3.13p we apply the standard p- approximation result in ID (see, 
e.g., [U Lemma 3.2]): there exists 4>p G 'Pp(^) such that 

W^U - MhH£) < Cp-^/^ ||^MIh3/2(^). (3.14) 

Putting together inequalities (|3.12|) - (j3.14|) and recalling our notation for the function (j), we 
obtain 

Wi^le - i^2,e\\o/ < Cp"^/^ IIV'I^ - V'2/|Ihi/2(£)- (3.15) 
Since ^ 

(u - uP) ■ h\QK = (u - ui - ) • h\QK = X] ^2/) e H^^^'^idK) 

with H^^^^idK) := {u e H-^/^{dK); {u, 1) o,dK = 0}, and the tangential derivative defines an 
isomorphism ^ : H'^/'^{dK)/R H~'^^^{dK) (see [HI Lemma 2]), we prove that 

Mi - V'2/|Ihi/2(^) <C ip- ^2/ 1/2, - <C\\{u- uP) ■ A||j^-i/2(a^). (3.16) 

ecdK ^ ' 

The first inequality in (j3.8p then immediately follows from estimates ()3.10p . ()3.15p . and (j3.16p . 
The second inequality in (j3.8p is true due to the continuity of the normal trace operator v — > 
V • h\QK as a mapping H(div, i(') H~^/'^{dK). □ 

Now, let us consider our Lipschitz surface T discretised by the quasi-uniform mesh. Using 
the Piola transform A4j for each element Tj and applying the ^-interpolation operator lip" on 
the reference elements, one can define the "global" H(div)-conforming /ip-interpolation operator 
Uf^ : H'-_(r) n H(divr,r) ^ Xhp (r > 0) such that for u'^p := Uf^u there holds 

Since the projection-based interpolation operator lip" preserves polynomial vector fields and 
provides H(div)-conforming approximations, the ^-interpolation error estimate (|3.7p on the 
reference element extends to the corresponding /ip-estimate in a standard way by using the 
Bramble-Hilbert argument and scaling. This result is formulated in the following theorem. 
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Theorem 3.1 Let u G H''_(divr, T), r > 0. Then there exists a positive constant C independent 
of h, p, and u such that 

||u - n^;,-u||H(divr,r) < Ch^'-^^'P^p-^ l|u||H'-_(divr,r)- 

4 Approximating properties of X/j^ in X 

The main purpose of this section is to prove that the orthogonal projection Php ■ X 
with respect to the norm in X satisfies an optimal error estimate (when the error is measured 
in the norm of X). The standard approach to tackle such kind of problems is to use the duality 
argument. Following this approach on F in our setting would require that for s G (0, ^] the space 
H~*(divr,F) (or, H~*(divr,F) if F is an open surface) is the dual space of H'*(divr,F) with 
respect to the H(divr, F)-inner product. However, this fact is true only for smooth surfaces. On 
the other hand, the duality argument can be applied face by face. This idea was suggested by 
Buffa & Christiansen in [11] and was successfully exploited by these authors in the context of 
the /i-version of the BEM for the EFIE. We will demonstrate below that, when employing the 
projection based H(div)-conforming interpolation operator, this approach provides an optimal 
/ip-approximation result in the energy norm of the EFIE. 

Let F be a single face of F. To simplify the presentation (in particular, to avoid imposing 
boundary conditions on the edges of dV which also belong to dV in the case that dT 7^ 0), we 
will assume that F is a (closed) Lipschitz polyhedral surface. All arguments extend to the case 
of an open piecewise plane Lipschitz surface in a straightforward way (cf. [llj). For the sake of 
simplicity of notation we will also omit the subscript F for differential operators over this face, 
e.g., we will write div for divj=,. 

Let H^(div,f) := {u e H^(f); divu G H''{r)} for r > 0. We wih also denote by Xftp(f) the 
restriction of X/^p to F. Then, for r > we introduce the operator Q^p : H''(div, F) X.fipit) 
as follows: given u G H''(div, F), we define QhpU £ X.hpiT) such that 

(u- QhpU, v)jj(^j^ p^ =0 Vv G X/jp(f) nHo(div,f), 

(4.1) 

QhpU ■ n = n^p^u • n on dT. 

Here, (•, •)H(div f) denotes the H(div, F)-inner product, and n is the unit outward normal vector 
to dT. It follows immediately from (j4.ip that 

||U - QhpU||H(div,f) < II" - nt"u|lH{div,f)- (4-2) 

Lemma 4.1 For r > /ei u G H'"(div, F). Then there holds 

1 /2 

||u - Q/.pu||H_v2(di.,f) < C ||u - n^>||H(div,f)> (4-3) 

where the constant C is independent of h and p. 
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Proof. We follow the technique used by Buffa and Christiansen in the proof of Proposition 4.6 
in [11] but rely on the properties of the H(div)-conforming projection based interpolation op- 
erator n^JJ'. For given r > and u G H^'(div, F) let us consider the following problem: find 
uq G H(div, r) such that 

(u-uo,v)jj(j;^P) = Vv G Ho(div,f), 
Uq • ii = n^J^u • n on df. 
Then, as an immediate consequence of Lemma 4.8 in [llj, there holds 

||u-uo||Hs+i/2(di^^f) <C||(u-n^;;u)-fi||j:^,(gp), s = -l, -i. (4.4) 
By the triangle inequality one has 

||u - QhpU||jj_i/2(ji^^p) < ||u - Uo||H-l/2(div,f) + ll"0 - QhpU|lH-l/2(div,f)- (4-5) 

Let us estimate each term on the right-hand side of (j4.5p . For the first term we have by using 
(fOl) with s = -I 

l|u-uo||H-i/2(div,f) < C'll(u-n^pu) -fill^.i^gp). (4.6) 

Denote u''^' := H^J^u. Since /^^(u — u^^) • ndcr = for any mesh edge £h C dV, we can apply 
the localisation result of Lemma 13.11 to estimate 

Let us fix an edge C dV. Then there exists an element K/^ C F such that = dV n dK^. 
This element Kh is the image of the reference element K under an affine mapping M, and let 
M be the corresponding Piola transform. Then using (|3.ip . (j3.2p and the standard property of 
the Piola transform (see [101 Lemma 1.5]), we have 

ii. hp^ ~u ((u-u^P) •n,(/;)o,£^ 

IKu - U '^) • n||^-l.^ = sup r—r 



A-1/2II-II 

= c/ii/2||(u-nfu).fi||^_,(,-), 

where i = M~^{ifi) C dK, u = MI~^{u\k^), if = ipoM, and n denotes the unit outward normal 
vector to dK. Hence, applying Lemma 13.31 and using standard scaling properties of the Piola 
transform (see [101 Lemma 1.6]), we estimate 



u — u 



• ^ C {if' II" - nf u||H(div,/.) < C [if' ||u - u''n|H(div,i.,)- (4-8) 
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Combining inequalities (j4.8p over all mesh edges £h C dF and recalling the notation for u''^, we 
deduce from (j4.6p and (|4.7p : 



1/2 ,. 

— div. 



11" - "o|lH-i/2(div,f) < C j l|u - nCu|lH(div,f)- (4-9) 
Now we focus on the second term on the right-hand side of (14. 5p and prove that 



||uo - Q.pu||H-i/2(di.,f) < C [1)'^' ||u - Uf;u\\^^,,^^^y (4.10) 

Denote Xp := Hg ^/^(div,f) and Ho/^(div,f) := Hi/2(div, f) n Ho(div, f). Let be the dual 

space of Xp. Then, by Lemma 4.7 in [llj, the operator / — V(div) : HQ''^(div, f ) — > Xj^ is an 

isomorphism. Moreover, it is easy to see that Qhp^o = QhpU and uq — QhpUo G Ho(div, F) C Xp. 
Therefore, we can estimate as follows: 

||U0 - QhpU||jj_i/2(di^_p) = ||uo - Q/ipUo||H-i/2(div^f) 

(v,Uo - QfepUo)x'-,X 



< C sup 



O^^vGX^ I|V||X^ 

(w - V(divw),uo - Q/ipUo)x'-,Xp 
< C sup ^ ^ 

0^weHj/2(div,f) l|w|lHl/2(div,f) 
(uo-Q/,pUo,w)h(j P) 

= C sup II — II ^ 

07^weHj/2(div,f) l|w|lHi/2(div,f) 

(UO - QftpUo, W - ndj;w)H(di,,f ) 

= C sup 



0^weHj/2(div,f) l|W|lHi/2(div,f) 

for the last step we used the definition of Qhp (see (|4.ip with v = II^p^w). Hence, using the 
Cauchy-Schwarz inequality and the interpolation error estimate of Theorem 13.11 (with L = F 
and r = ^) we prove that 



||uo - Qhp^n-y ^idi.,r) ^'^[pj ll"o - 2/xpu|lH(div,f)- (4-11) 

The norm on the right-hand side of ()4.1ip is estimated by applying the triangle inequality, ()4.4p 
with s = — |, (j4.2p . and the continuity property of the normal trace operator: 

||uo - Qhp^n{div,r) ^ ll"o - u||H(div,f) + II" " Qhpu\\u{div,r) 

< c\\{u- uf;u) ■ n||^-i/2(af) + 11" - n^>llH(div,f) 

< C||u-n5i;;u||H(di.,f)- (4.12) 
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The desired inequality in (j4.10p then follows from ()4.1ip and (j4.12p . 

To obtain (j4.3p it remains to collect (j4.9p and (j4.10p in (j4.5p . This finishes the proof. □ 

Now we can prove the main result of this section. 

Theorem 4.1 Let Php '■ X X^p be the orthogonal projection with respect to the norm in X. 
If u £ XJ' with r > — ^ , then 

with a positive constant C independent of h, p, and u. 

Proof. First, let us assume that r > 0. We consider the discrete vector field v G X/^p such 

that v|p(i) = Q^p(u|p(i)) for each face T^*) of T (here, Q^^ denotes the operator defined as in 

(j4.ip with respect to the face F*^*)). Due to the localisation properties of the norms in X and in 
H(divr,F), we have by Lemma |4. II 

I 

||u-Pftpu||x < ||u- v||x < C ^||u|rw - Qip(u|rw)llH-i/2(div^(^),rW) 

i=l 

- ^ (p) II" ~ "llH(divr,r), 

and inequality (j4.13p then follows from the error estimate of Theorem 13.11 

Now, let r G (— ^,0]. Assume that u G X** with some s G (0,so)! where sq G (0, |] is the 
same as in Lemma 13.21 Then, using the first part of the proof, one has 

||u-P,,pu||x < C(Jj ||u||x- 
On the other hand, it is trivial that 

||U - PhpU||x < ||u||x- 
Therefore, applying the interpolation argument which relies on Lemma 13.21 we prove 

||u-P,.pu||x<C(|j ||u||x^ VUGX^ 

This estimate yields (|4.13p due to the density of regular functions in X**, and the proof is finished. 

□ 
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